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Abstract 

Renormahzed Hamiltonians for gluons are constructed using a perturbative boost-invariant 
renormahzation group procedure for effective particles in light-front QCD, including terms up 
to third order. The effective gluons and their Hamiltonians depend on the renormahzation 
group parameter A, which defines the width of momentum space form factors that appear 
in the renormahzed Hamiltonian vertices. Third-order corrections to the three-gluon vertex 
exhibit asymptotic freedom, but the rate of change of the vertex with A depends in a finite 
way on regularization of small-x singularities. This dependence is shown in some examples, 
and a class of regularizations with two distinct scales in x is found to lead to the Hamiltonian 
running coupling constant whose dependence on A matches the known perturbative result from 
Lagrangian calculus for the dependence of gluon three-point Green's function on the running 
momentum scale at large scales. In the Fock space basis of effective gluons with small A, the 
vertex form factors suppress interactions with large kinetic energy changes and thus remove 
direct couplings of low energy constituents to high energy components in the effective bound 
state dynamics. This structure is reminiscent of parton and constituent models of hadrons. 
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1. Introduction 

Current studies of liadronic structure are guided by three physical pictures. The first picture 
is based on the constituent quark model, which serves classification of hadrons in particle data 
tables [|[]. Quantum numbers of a hadron in this model correspond to a simple Hamiltonian 
with only kinetic energy of two or three quarks and inter-quark potentials in the hadron rest 
frame, with no gluons. The second picture is provided by the parton model for hadrons in 
the infinite momentum frame P]. Modern versions of the model introduce a slew of quarks 
and gluons with distribution functions in variable x - a fraction of the hadron momentum 
that is carried by a parton. About half of the hadron momentum is carried by gluons, with 
mostly small values of x, so that many partons can share the hadron momentum. Binding of 
partons is not described by the parton model. In the third picture, hadrons are considered to 
be excitations of a complicated ground state (vacuum) that contains condensates of quarks and 
gluons. Understanding of hadronic structure in the third way relies on the assumed ground 
state properties [Q. Despite recent progress in experimental and theoretical studies of hadronic 
structure, including lattice approach 0, |^, the three basic pictures are not yet unified in a 
single quantitative formulation of QCD. To connect constituent quarks and partons with QCD 
degrees of freedom, one needs a relativistic description of effective particles in quantum field 
theory. 

This paper describes a perturbative third-order calculation of renormalized Hamiltonians 
for effective gluons in the light-front Fock space. The effective gluons are derived in a boost- 
invariant renormalization group procedure for particles p, which originates in the similarity 
approach to renormalization of Hamiltonians [0 and the notion of vertex form factors for ex- 
tended strongly interacting particles |^. The renormalization procedure provides a connection 
between the canonical quantum field theory and the concepts of bound state constituents in 
the rest and infinite momentum frames. For simplicity, this paper is limited to gluons. Quark 
effects in the gluon dynamics are mentioned only in passing. Gluons alone are worth a dis- 
cussion since their interactions are responsible for asymptotic freedom. This feature requires 
understanding in Hamiltonians independently of quark dynamics. Also, asymptotically free 
effective gluon interactions display specific sensitivity to regularization of small-x singularities. 

Section 2 presents the initial regularized Hamiltonian for gluons. The Hamiltonian includes 
ultraviolet counterterms that are calculable order by order in the procedure described in Section 
3. The procedure introduces vertex form factors in the effective gluon interactions. The form 
factor width parameter. A, is reduced from infinity down to the scale of hadronic masses through 
a solution of a differential equation, which eliminates large momentum transfers from the bound 
state eigenvalue problem for effective gluons with small A. The coupling strength of the three- 
gluon vertex as function of A, is calculated in Section 4 and analyzed in Section 5. These 
two Sections show how asymptotic freedom of effective gluons emerges in the light-front Fock 
space Hamiltonians for QCD. Section 6 provides a short summary and a brief discussion of 
how the effective particle calculus can be applied to electron-hadron scattering in a simplest 
approximation. 

2. Initial Hamiltonian 

The canonical light-front QCD Hamiltonian requires regularization and counterterms 0. 
To regulate the Hamiltonian, momenta = + and p^ = {p^-,p^) are parametrized using 
the + momentum ratios x and relative transverse momenta that will be described below. 
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Regularization is imposed through factors that exclude large \k^\ and small x, preserving all 
kinematical symmetries of the light-front dynamics (i.e. the Poincare symmetries of the surface 
= x^ + x^ = in space-time), and processes of creation of particles from the bare vacuum are 
absent. Power-counting and renormalization strategy for the absolute co-ordinates, x~ = x^—x^ 
and x"*" = (x^, x"^), or p'^ and p"*" [^, are modified when one goes over to the variables x and k"*". 
However, key features remain similar and perturbative results described in the next Sections 
agree with expectation that the ultraviolet renormalization of light-front Hamiltonians involves 
functions of x. 



2. a Canonical terms 

The classical Lagrangian density for gluon fields is 

C = -hrF^^F"^ , (2.1) 

where F^^ = d'^A" - d^A'' + iglAi", A"], and = A'^H'' with [f, t''] = if'H". The Lagrangian 
implies equations of motion, d^F^^^ = ig[F^'^ , A^j], and for fields satisfying these equations the 
canonical energy-momentum density tensor is T'^'^ = —jr<^t^'^g'^j['^ -)_ gt^'^F"''^^F^p/4:. 

In the gauge A""^ = 0, the Lagrange equations constrain A" to A" = {l/d^)2d^A-^ — 
(2/d'^'^)ig[d~^A^, A^] and the independent field degrees of freedom are A-^. The first term 
in A~ is independent of the coupling constant g, and can by definition be included in a new 
constrained field A^^ = [A^ = 0,A~ = {l/d^)2d^A-^, A-^]. The second term can be kept 
explicitly as part of interactions. Using this convention and freely integrating by parts, one 
obtains an expression for light-front energy of the constrained gluon field 

P- = \ [ dx-d^x^n\^+=o , (2.2) 
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where Ti. = T+ and 



with 



+ + Ha* + 'HyQAA]^ , (2.3) 



nA^ = -lA^{d^YA^ , (2.3a) 

nA^=gido.A''^[A",A''r , (2.36) 

HA^ = -^g'[A^,ApnA'^,Af'r , (2.3c) 

n^aAAP = [id'-A^, A^r-^^[id+A^, A'-r . (2.3rf) 



This expression is a candidate for further consideration in analogy to QED |jTO|, |TT|, |T^, |T^. A 



heuristic expression for the quantum gluon energy operator is obtained by substitution 

^' = Y.j \k] [t^sla,^,e-^'^ + feZal^e^'^] ^^^^ , (2.4) 
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where [k] = 9{k~^)k~^<Pk-^ / (167r^A;+) and e^.^ = {e^^ = 0, = 2k-^e-^/k~^, e-^). a numbers gluon 
spin polarization states and c is a color index. The creation and annihilation operators satisfy 
commutation relations 

= k+ ^{k - k') 5""' 5"'' , (2.4a) 

where ^{p) = 16TT^6{p'^)6{p^)6{p'^), and commutators among all as, and among all a''"s, vanish. 
For all momenta, spins and colors, akad^) = and al^^ creates bare gluons from the state |0). 

The plain insertion of Eq. (2.4) into produces terms with creation and annihilation 
operators appearing in all possible orders. All terms are then ordered so that creation operators 
stand to the left of annihilation operators. The process of ordering produces commutators of 
creation and annihilation operators, which lead to diverging integrals. All such terms can be 
dropped at this stage entirely since they will be either removed by regularization, in the case 
of modes with k~^ = or, after regularization, they will be replaced by well defined mass 
counterterms that result from a renormalization group procedure and contain free additive 
finite parts. 

The ordered operator, denoted Pquantum^ is highly divergent. For example, a correction order 
g"^ to the free energy k~ = k-^^/k^ of a single bare gluon state \kac) = al^^O), diverges due 
to integration over an infinite range of transverse momenta of virtual gluons that appear in 
the intermediate states of second order perturbation theory. The energy correction diverges 
also due to small x singularities. Namely, the gluon momentum k~^ can be shared by two 
intermediate gluons carrying fractions x and 1 — x. The sum over intermediate states involves 
an integral over x from to 1 while the polarization vectors of intermediate gluons provide 
x and 1 — a; in denominator of the integrand, cf. in Eq. (2.4). As another example, the 
product PquantumPquantum IS evcu morc divergent than the energy correction because it does 
not contain the energy denominator that reduces the contribution of intermediate states with 
large momenta in perturbation theory. Consequently, exp {—iPquantum^^ 1'^) ^ candidate for 
a unitary evolution operator in time x"*", is not defined before one regulates Pquantum by limiting 
the range of momentum that the bare gluons may have. 



2.b Regularization 

The first step in regularization procedure is made by limiting the range of momentum 
integration in Eq. (2.4), cf. Ref. 0. Let \k-^\ < Vt and /c"*" > e"*", with understanding that 
Vt ^ oo and e+ ^ when the regularization is being removed. 

The lower bound on k^ implies that the regulated expression for Pquantum^ denoted by Pq^+, 
does not contain any terms with exclusively creation, or annihilation operators. Such terms 
would be forced by a translationally invariant integral over x~ to preserve momentum A;"*", 
while the momentum they would have to create, or destroy, is at least ne"*", where n denotes the 
number of creation, or annihilation operators in such terms, respectively. These two conditions 
are incompatible. Hence, P^^+ does not contain terms that could alter the bare vacuum state 
|0). The coupling constant g is assumed to be sufficiently small for stability of the regularized 



theory built on top of |0), cf. ||14 



The limits on absolute momenta in -P^eH- violate boost invariance of light-front dynamics, 
and in the term (2.3d) with inverse powers of 9^, one may still have in denominator. To 
eliminate the violation of boost invariance and regulate the singularities, P^^+ is further 
curbed through the following step. 
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Interaction terms in Pq^+ are modified so that changes of the particle momenta are limited. 
In this work, the transverse momentum changes are limited by a parameter A, and changes of 
X by a parameter S. As an example, it is useful to consider the three-gluon vertex, 

HAsne+ = E / [123] m + k2- h) [Cn.+ (123) alala^ + C*,+ (123) ala^a^] , (2.5) 

123 

Momentum conservation implies that one can write = Xik^, = Xik^ + k^2j — ^2^3"? 

= X2k^ — Hi2y with X1+X2 = 1, and = X2ki — Xifc^. Xi, X2, and ^,^2 invariant under 
seven kinematical transformations of light-front dynamics. It is helpful to call the momentum 
carried together by all annihilated or created particles in a single vertex a parent momentum in 
the vertex. In the vertex (2.5), k^ is a parent momentum. Also, a slash in a subscript is used 
below to indicate that the momentum before the slash sign is considered to be a daughter of the 
parent momentum after the slash sign. For example, ki = xi/^k^ + /t^yg, where xi/3 = xi/x^ 

and ^]^^3 ~ ^12* 

The momentum changes in the vertex (2.5), are limited by inserting a factor r^sii^j/i^Xi/s) 
for each creation and annihilation operator. Since k^ is a parent momentum, 0:3/3 = 1, and 
^3/3 ~ 0; regularization factor for the term (2.5) equals rA<5(K"'"^, x)rA<5(K-'-^, 1 — x), 

where x = xi/3 and k"*" = /t^2- Factors r^s are chosen to have the form 

TAsi^^^x) =rA{t^^y5{x)9{x) , (2.6) 

where 

rA(^) =exp(-^/A2) , (2.7) 

and rs{x) suppresses the region of x smaller than 6. That r^lz) falls off exponentially is 
a guarantee for ultraviolet convergence of all transverse momentum integrals that appear in 
perturbation theory. Integrals that behave as In A or A*^ with positive n for A 00, will 
be called ultraviolet divergent. The small-x regulating function, rs{x), must vanish sufficiently 
quickly for x — to regulate all small-x singularities. In a sense to become clear in Section 5, 
the factors rs considered in this work lie in the vicinity of two cases, 

rsix) = x/{x + 5) , (2.8) 

and 

rs{x) = e{x-6) . (2.9) 

Integrals that behave as In 6 or 5" with negative n, will be called small-x divergent. Mixing of 
the ultraviolet and small-x regularizations through expressions of the type exp [— (k-'-^/x)/A^], 
is discussed in Appendix E. 

Every creation and annihilation operator in all vertices in the operator Pq^+ is supplied with 
a factor r^s, k^ being replaced by the corresponding parent momentum. Terms that contain 
four operators are recast as contracted products of terms with only three operators, which are 
already regulated. Details of the quartic terms are given in Appendix A. The instantaneous 
terms containing inverse powers of are regulated in the same way, by interpreting the 
momentum that is transfered along the inverted as a momentum carried by a virtual particle 



that connects two vertices ||Tl|, |T^. The fully regulated operator P^^+ is denoted by 



A<5 
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In the last step of defining the initial Hamiltonian, H^s, the limits of Q 
are taken with A and 6 kept constant, 



oo and e'^ 



AS 



lim lim 



AS 



AS 



(2.10) 



X/^s denotes counterterms, which need to be found. The renormalization group procedure that 
provides means for finding them order by order in perturbation theory, is described in next 
Sections. The initial Hamiltonian has then the form, cf. Eq. (2.3), 



Has = Hji2 + Hj^a + Hj^i + H^qj^a]^ + X^s 



where, for example. 



and 



Has = J2 [123] ^{p'^ - p) rA5(3, 1) \g F123 aja^ag + g Y*^;^ 4a2ai 



123 



(2.11) 



(2.12a) 



(2.126) 



These two terms are quoted from the Appendix A, where all regulated canonical terms are listed 
and the notation is explained. X^s will be discussed below and in next Sections. Note that the 
free Hamiltonian (2.12a) contains no regularization. This is necessary to preserve kinematical 
light-front symmetries. Differences between the bare three-gluon vertex (2.12b) and effective 
gluon vertex from Section 5, are described there. 



3. Renormalization group procedure for particles 

In Eq. (2.11), H^s is expressed in terms of the creation and annihilation operators for bare 
gluons in regulated local theory. This Section describes the renormalization group procedure 
that is used in Section 4 to re- write the initial Hamiltonian H^s in terms of operators that create 
or annihilate effective gluons, instead of bare ones. The effective gluon operators are obtained 
by applying a unitary transformation L(\ to the initial bare operators. The effective operators 
depend on the parameter A that labels Ux. 

The parameter A has dimension of mass and distinguishes different kinds of effective gluons 
according to the following rule. Effective gluons of type A can change their relative motion 
kinetic energy through a single effective interaction by no more than about A. The transforma- 
tion I4x is mathematically designed in perturbation theory so that resulting interaction terms 
contain vertex form factors and the latter limit the kinetic energy changes by their width pa- 
rameter, which equals A. All Hamiltonians with different As are equal and the re- writing does 
not introduce any change in the theory, although the same Hamiltonian appears differently 
when expressed in terms of different gluons. For brevity, the effective gluons corresponding to 
some value of A are referred to as gluons of width A. 



3. a General features 

The form factor width parameter A greatly differs from the regularization cutoffs, because it 
may be kept finite, even small, while the cutoffs have to be made extremely large to approximate 
the initial theory. Even if the expansion in terms of bare particles is hopelessly complicated, a 
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hadron may still have a well defined, convergent expansion in the basis of effective constituents 
with small width A. 

Thanks to the vertex form factors, in the Fock space basis built from gluons of width A, the 
effective Hamiltonian matrix elements quickly tend to zero when the effective gluons change 
kinetic energy across the matrix element by more than A. Therefore, the effective Hamiltonian 
matrix is narrow. This is important for application to bound state physics because eigenstates 
of narrow matrices may have a small number of dominant components [|15|. In the case of 
hadrons, the constituent model suggests that the intricate complexity of QCD is buried in the 
structure of constituents and their interactions, while the number of effective constituents is 
small. The success of perturbative QCD in reproducing changes of deep inelastic structure 
functions with momentum transfer down to fairly small values, suggests that the structure of 
effective constituents can be approximated using perturbation theory. The idea of effective 
particles is by no means new [|l^]. New element is the renormalization group procedure for 
effective particles in QCD. 

The renormalization group procedure provides means to find counterterms X^s that have to 
be included in the initial Hamiltonian to compensate for the spurious effects of the ultraviolet 
regularization. One takes advantage of the form factors in the Hamiltonian vertices in analogy 
to Ref. 0. The narrow dynamics can smear states of effective particles with finite energy, only 
by less than nA towards high energies in nth order perturbation theory. For to raise the free 
energy of a state by nX, the effective interaction must act on the state about n times when the 
form factors die out exponentially for energy changing by more than A. The highest order n 
that is still independent of A, approaches oo when A oo. Therefore, to obtain the ultraviolet 
regularization independent results, at least in perturbation theory to all orders, it is sufficient to 
demand that the Hamiltonian coefficients in front of creation and annihilation operators with 
finite A are independent of the ultraviolet regularization. Hence, one can read the structure of 
ultraviolet counterterms from the coefficients, see Eq. (3.11). However, the effects of small-x 
regularization are not under control of the renormalization group procedure and the coefficients 
with finite As may depend on rs{x). The dependence on 5 — sometimes drops out but finite 
effects may remain, as will be shown on examples in Section 5. 

In the perturbative renormalization group procedure for deriving effective particles and 
their interactions, one never encounters genuine infrared singularities associated with small 
energy denominators. This is explained below Eq. (3.10), where the differences of invariant 
masses are taken care of in analogy to differences of energies in the similarity renormalization 
group procedure for Hamiltonians [0]. The perturbative denominators are effectively limited 
from below by A. The non-perturbative part of dynamics with relative motion kinetic energy 
changes smaller than A, is first tackled when one proceeds to solve the effective Schrodinger 
equation. Since the form factors keep the effective dynamics in a well defined range of energies. 



numerical methods may apply in finding approximate solutions to the full theory |15 



3.b Construction of H\ P 

Let a commonly denote the bare operators, akac or a|,^^. Operators a are transformed by 
the unitary operator h(\ into operators a\ that create or annihilate effective particles of width 
A, with identical quantum numbers. 

ax = UxaUi. (3.1) 
The initial Hamiltonian H^s is re- written in terms of a\, H/^s = Hx{a\)- If quarks were 
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included, i^AS would correspond to the QCD Hamiltonian written in terms of canonical quarks 
and gluons, associated with bare partons, or bare currents. H\ for A comparable with masses 
would represent the same Hamiltonian written in terms of constituent quarks and gluons. 
Applying Ux, one obtains 

nx = Hx{a)^UiHA5Ux . (3.2) 

Hx has the same coefficient functions in front of products of as as the effective Hx has in front 
of the unitarily equivalent products of a^s. Differentiating Tlx with respect to A, one obtains 

n', = -[rx,nx] , (3.3) 

where Tx=U'lUx. 7a is constructed below using the notion of vertex form factors. For example, 
if an operator without form factors has the structure 

dx = j [kik^h] Vx{l, 2, 3) a\^^a\^,axk, , (3.4) 
the operator with form factors is written as fxOx and has the structure 

fxOx = I [hk2h] fx{Mi2, Ms) Vx{l, 2, 3) ai,^ai,^axks , (3-5) 

where 

fxiM^2, Ms) = exp[-(M?, - Ml)'/X')] . (3.6) 

For any operator O expressible as a linear combination of products of creation and annihilation 
operators, fO contains a form factor fx{Mc: Ma) in front of each product, where Mc and Ma 
stand for the total free invariant masses of particles created (c) and annihilated (a) through 
the product, respectively. For gluons in Eq. (3.6), Ms = 0. 

The relative motion kinetic energy changes in interaction vertices of effective particles, are 
limited by demanding that Hx = fxGx, with some unknown Gx- Then, one derives equations 
for Gx that result from the choice for fx and some definition of 7^. In practice, one first uses 
Ux to transform the Hamiltonian to TCx = fxQ\ and then one calculates coefficients of as in Qxi 
which are the same as the coefficients of ax in Gx- This calculation includes the construction 
of Txi and proceeds as follows (the subscript A is dropped for simplicity of notation). 

n' = fg + fg' = -[r, g^] - [r, g:\ . (3.7) 

g is split into two parts. The free part is bilinear in as and independent of the coupling 
constant g. The remaining, interaction dependent part is denoted by In the present work, 
^0 is taken to be independent of A. Eq. (3.7) contains two unknowns, T and Without loss 
of generality, one assumes that T ^ when t// ^ 0, and one expands operators in powers of 
gj with the goal of enabling the procedure to work order by order. The expansion of T starts 
from the term of order g^. Changes of ^/ with A should start from second power. If 

fg'r = -f[T, gi] , (3.8) 

then 

[T, go] = [(1 - f)gi]' , (3.9) 
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and 



G', = [fGj, {(1 - /)^,}go] ' (3-10) 

where the curly bracket with subscript Qo, indicates the solution for T that follows from Eq. 
(3.9). The choice of / made above implies that perturbation theory for T and Q does not lead 
to small energy differences in denominators, since 1 — / vanishes quadratically with the energy 
difference. Qj contains only connected interactions because Eq. (3.10) has a commutator on 
the right-hand side. The initial condition for Eq. (3.10) is provided by H^^s, so that Eq. (3.10) 
in the integral form reads 



Has + 



ds 



fsGis, {(1 - fs)Gls} 



Go 



(3.11) 



which allows one to find the counterterms X^s using the condition that they remove dependence 
on regularization from the second term for finite A and relative momenta of interacting particles. 
Finally, = hOx- 



3.C Perturbation theory for Qjx 

This Section contains formulas used in Section 4 for solving Eq. (3.11) in perturbation 
theory up to third order. Since the perturbative expansion has formally the same structure as 
in scalar theory [|l^], only main steps are listed. In the first step, Qix is expanded into a series 
of terms Tn ^ g"", 

oo 

Gl = J2^n . (3.12) 

n=l 

It immediately follows from Eq. (3.11) that ri is independent of A and equal to the second 
term in the initial Hamiltonian from Eq. (2.11), i.e. (2.12b). One has ti = a2i + ai2, where 
a2i denotes the first and au the second term on the right-hand side of Eq. (2.12b). The left 
subscript denotes the number of creation and the right subscript the number of annihilation 
operators. The corresponding Hamiltonian interaction term is obtained by multiplying the 
integrand in Eq. (2.12b) by fx = exp [—{ki + k2)^/X'^] and transforming as into axs. 

For T2 = Pn + p3i + As + /322, Eq. (3.10) implies 

^2 = [{Ai},/ri] = /2[riri] , (3.13) 

where /2 = {/'}/ — /{/'}? with understanding that the first factor / in all terms of /2 is for 
the first r in the square bracket, and the second factor / in all terms of /2 is for the second r 
in the bracket. The square bracket denotes all connected terms that result from contractions 
replacing products Oja^ by commutators [aj,a]]. The solution for T2 is then 

T2A = ^2x[riTi] + r2oo , (3.14) 

where ^2A = fi depends on incoming and outgoing momenta in the two vertices formed by 
the operators in the square bracket. In the sequence aTabbrbcC, the three successive configurations 
of particle momenta are labeled by a, b and c. To write down compact expressions for J^2X, 
the symbol uv = J^'^v ~ -^lu defined, where A^^„ denotes the free invariant mass of a set 
of particles from configuration u that are connected to the particles in configuration v by an 
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interaction r^^ in the sequence uTuvV. Spectators of the interaction r^^ do not count. In this 
notation, 

fx{Mab,Mba) = exp [-{abyX^)] = fat • (3.15) 

The parent momentum for the vertex connecting two configurations u and v, is denoted by 
Puv, and in the fohowing equations Puv is written in place of P^. In all expressions, the minus 
component of momentum of every gluon is given by the eigenvalue of Qq — Hj^2 from Eq. 
(2.12a), i.e. k~ — k^'^/k'^. Thus, with the Gaussian vertex form factors, 

•^2(0, b, c) = + fcc^ if-bfbc - 1] ■ (3.16) 

In Eq. (3.14), 

where Xas2 denotes all ultraviolet counterterms proportional to g^. 

For third order terms, T3 = 721 + 712 + 741 + 7i4 + 732 + 723, Eq. (3.10) gives 

^3 = [M, {(1 - /)r2}'] + [{/Vi}, fr,] . (3.18) 

After integration, 

T3A = ■^31a[ti[tiTi]] - .F32A[[Tiri]Ti] + ^2a[t2ooTi + riT2oo] + T300 , (3.19) 

where, for any sequence aTabbrbcCTcdd, 

^32(0, 6, c, d) = -^3i(d, c, 6, a) , (3.20) 



J^3i{a,b,c,d) = 



Pcbcb + pcdcd 
cb^ + cd'^ 



iPbdbd + Pbaba) 



fabfbcfcdfbd 1 



fabfl 



bd 



ab^ + bc^ + cd^ + ah^ + bd'^ 



+ 



+ P6d 



bc^ + cd'^ 



db 



bd 



1 



fabfbcfcd 1 

a62 + + cd^ ' ah^ + bc^ + cd^ + bd^ 



fabfbcfcdfl 



(3.21) 

The last term in Eq. (3.19) denotes counterterms proportional to g'^, i.e. rsoo = Xass- The 
next Section describes calculation of 721. The calculation requires knowledge of P22 and 
/Ssi. For all terms, TTij = 7r]j. 



4. Interactions of effective gluons 

This Section describes derivation of effective gluon dynamics in third order perturbation 
theory in the coupling constant gx that measures the strength of effective three-gluon vertex 
in the Hamiltonian Hx{ax)- Since g\ = g + o{g^), Hamiltonian terms of order gx and gf^ can 
be calculated using expansion in powers of g before one proceeds to the third order terms that 
define gx- All calculations are carried out in the framework described in Section 3. 

In the re- written Hamiltonian Hx{ax), the terms that have coefficients of order 1, are 



10 



H{0) =Y. W\^Akac<^^kac ■ (4.1) 

The subscript A indicates that a\kac annihilate and a\j^^^ create effective gluons of width A. 
The effective gluons can also be interpreted as having a spatial transverse width on the order 
of 1/A for moderate values of x. This interpretation is explained below Eq. (4.3). 

The terms in H\{a\) that have coefficients order are (see Appendix A for details of the 
notation) 



123 



123 O'XzO'XlO'Xl 



(4.2) 



where 



rs{x) =rs{x)rs{l- x) , (4.3) 

and the form factor /a(-Mi2,0) = exp [— K;j'-2*^/(a;ia;2A^)^] falls off as function of the relative 
transverse momentum at a rate that depends on x carried by gluons. For moderate values 
around 1/2, the transverse momentum width is on the order of A/2, but for x approaching 0, 
the transverse momentum width of the vertex becomes very small, leading to a spread of the 
interaction strength in the transverse spatial directions. Thus, the coupling of effective gluons 
to the wee region is quite different from the canonical coupling in Eq. (2.12a). 

Terms with coefficients order g"^ are derived by changing a to a\ in t\2 = Pxn + /3a3i + 
P\i3 + P\22, and by inserting form factors as described in Section 3.b. The contribution of T200 
includes a counterterm induced by the ultraviolet regularization r^. Namely, to evaluate terms 
Pxii one needs to know the counterterm jSoou- It follows from Eq. (3.14) that 

pMi = i:J[k]§al^a,.^ , (4.4) 

where 

lA = i^j'^^Jd'^^^:F2x{KK,k)2Y.\Y,2u\'fl, + , (4.5) 

lOTT'^ JO X(l — Xj J k+ ^ 

and the last term, /x^, is contributed by the counterterm /3ooii- The structure of /3ooii is known 
here from hindsight, i.e. the regularization dependence of the integral in Eq. (4.5) results in 
a number that depends on the function r/^{H?) but does not depend on the gluon quantum 
numbers. One has k~ = k^'^/k'^, = (TW^ + k^'^)/k^, = k-^^/[x(1 — x)] and rk,i is 
given at the end of Appendix A. The sum over quantum numbers of intermediate two-gluon 
states is E12 l>"i2fcP = NcK^[1 + 1/x^ + 1/(1 - xf] = k'^P{x)/[2x{1 - x)], where P{x) is the 
Altarelli-Parisi gluon splitting function Pgg{x) [|19|. A^^c = 3 denotes the number of colors. 

Assuming that for some A = Aq the effective gluon mass squared should have the value /Iq^, 
one can calculate the counterterm mass //^ from Eq. (4.5). The resulting effective mass takes 
the form 

= f^ls + -£, / / d'-^ ip[^2A - -F2aJ .'P{x)/[xil - X)] , (4.6) 

which is independent of the ultraviolet regularization. 



II 



rs,{x)^[h{x)f^T]{x)Tl{l-x) . (4.7) 

To find a more convenient notation for tlie arbitrary, possibly 5- dependent number //q^ 
can be replaced by an expression that follows from a second order perturbative result for a 
single effective gluon mass correction obtained from the eigenvalue equation for Hx{ax). This 
expression introduces a new parameter fij that would equal the second order result for physical 
gluon mass squared, if such gluons existed. With the form factor / of Eq. (3.15), one obtains 

/^A = + £ dx rs^,{x)P{x) dz exp [-2^7 A"] , (4.8) 

while the counterterm mass is 

IjI^ = + j^dxr5^{x)P{x) dze^v[-^zx{l-x)/A^] . (4.9) 

Inclusion of n/ flavors of quarks with bare masses m/, produces an additional term in the 
integrand in Eq. (4.8), equal Z]/=i Qf{^) ^xp [—2z'j/\'^] with Qf{x) = 1 — 2a;(l — x) + 2m'j/zf 
resembling the gluon splitting function into a pair of massive quarks of flavor /, and Zf = 
z + mj/[x{l — x)]. Quarks would also add Yl'fLiQfix) to P{x) in Eq. (4.9). Note that for 
convergence of the integral over x in the effective mass term (4.8), it is sufficient that the 
regulating function behaves as rs{x) ~ a;^ with e > 0. In the counterterm (4.9), rs{x) has to 
vanish at small x at least as fast 

Having the result for P\u, one replaces the bare a with effective ax. The second order 
effective gluon mass term is 

2 

^(2)ll=E/Wf^«Lc«Afc.c . (4.10) 

Other terms of order g"^ are derived following the same path, but they do not require ultraviolet 
counterterms for regularization adopted in Section 2. 

Of main interest to this work is the term ■yx2i whose calculation involves Ti, (322 and 
P31. Details of the calculation are described in Appendices B and C. Eq. (Bl) for ^X2i has the 
following structure, cf. Eq. (3.11), 

7A21 = 7oo2i + / dxd^ti^ [FxC + D]Ra . (4.11) 

Fx denotes form factors that multiply terms C, which means that momentum integrals in terms 
C cannot produce dependence on the ultraviolet regularization factors Ra in the limit A 00. 
The first step in calculating 7^21 is to evaluate terms / DR^. and find their dependence on 
Ra, to construct the counterterms 70021 ■ This step is carried out in Appendix C, including n/ 
flavors of quarks, with the result, see Eqs. (B9) and (C18), that 

7oo2i = E / [123] m + k2- h) 7oo a\alas , (4.12) 

123 -'-"^ 

where 

loo ^Yi23^^n-{N,[ll + h{xi)]-2nf} + -f finite , (4.13) 

O Li 
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and 



6 / dxrstix) 



1 



1 — X X — Xi 



1 



X 



9rs{xi) / dxrsf,{x) 



1 



1 



1 — X. 



(1 



(4.13a) 

f5(xi) is given by Eq. (4.3), rs^{x) by Eq. (4.7), and rst{x) by Eq. (B4a). The counterterm 
part finite removes finite effects of ultraviolet regularization. It contains unknown functions 
of Xi that multiply structures Y12, 1^13, Y23 in ^123, and the additional structure Y4, see Eqs. 
(CIO) and (C15). 'j finite differs from canonical terms but it does not influence the third order 
running of the Hamiltonian coupling constant. 

Having found the counterterm, one obtains a finite answer for 7^21 in the limit A — 00, 



7A21 = E / [123] #(p^ - p) [WX12Y12 + Wxi3Y,3 + WX23Y23 + WxaY,] al4a3 , (4.14) 

123 

where the four coefficients W\ are functions of Xi and ^12^- The complete expression for 7^21 
is given in Appendix B. The effective three-gluon interaction term of order in Hx{ax) is 
-f^(3)2i + -f^(3)i2, where -^(3)12 = -^^(3)215 and 

Hi3)2i = fxUxlx2iUi . (4.15) 

The operator lAx is unitary with accuracy to terms of order higher than 3rd and its action is 
equivalent here to changing as to aAS as in lower order terms. 



5. Three-gluon vertex 

This Section describes the running coupling constant gx that appears in Hamiltonian in- 
teraction vertices written in terms of gluons of width A. The word "running" means changing 
with A. The coupling constant is defined through the strength of the three-gluon vertex for 
some value of xi = xq when — > 0, in direct analogy to the coupling constant in 0^ theory 
in 6 dimensions ||18| or electric charge in the Thomson limit in QED. The effective three-gluon 
vertex is a sum of ^^(1+3)21 and -f^(\+3)2i) where i^(i+3) = H{i) + H(^') and 



^^(1+3)21 = E / [123] ^{P^ - P) fx [VXI2Y12 + VaisFis + yx23>^23 + WxaY,] a\A2<^X3 , (5.1) 
123 

with fx = exp — [ni2^/ (xiX2A^)]^. The effective vertex contains creation and annihilation opera- 
tors for effective gluons of width A, instead of bare ones in Eqs. (2.12b) and (4.14). This change 
is reflected by the presence of the form factor fx, which strongly suppresses emission of effec- 
tive gluons with transverse momentum larger than X1X2A. The vertex functions VAi2(a^i, '^12^)5 
^Ai3(2^ij ^12^) and VA23(a;i, '^12^) are used below to evaluate the running coupling constant. The 
fourth function, Wx4,{xi, 1^12^), multiplies I4 which is distinct from canonical structures. Wa4 is 
independent of A in the limit k]^2 ^ and does not contribute to the running coupling constant. 

When K]^2 ~^ 0) all three vertex functions Vxij for ij = 12, 13, 23, vary with A equally. Using 
results from Appendix D without quarks, one can introduce a single function, 

Wx{x) = Vx^j{x,0^)-Vx,^J{x,0^) 
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iVc [n + h{x)] In A , (5.2) 



487r2 " ^ ' ^ Ao 

to describe behavior of all three functions V\ij{x, 0^) = grs{x) + Wx{x) for finite x and A in the 
limit 5^0. 

W\{x) depends on x through h{x). The case h{x) =0 corresponds to the standard asymp- 



totic freedom result |]T6| in the following sense. At certain value of A = Aq, the three potentially 
different numbers gxoij = Vxoij{xo, 0^) can be set equal to a common value go = gxo by adjust- 
ing finite parts of the counterterm (4.13), to match the gauge symmetry result from the bare 
canonical Hamiltonian of QCD. For h{x) = 0, the choice of Xq does not matter as long as 6 
is negligible in comparison to Xq and 1 — Xq. Then, for A 7^ Aq the three couphng constants 
remain equal, gxij = Vxij{xo, 0^) = g^ = gQ + Wx{xo) + o{g^), with g = g^ in Wx{xo), or 

gx = go - -^llN.ln^ . (5.3) 

Expressing go by gx in the differential Eq. (3.10), or in its integrated forms, (3.11) or (5.3), in 
the perturbative expansion up to third order in gx, one has 

A^(7A = /3ogl , (5.4) 

with 

ft—i^. (5.5) 

which is equal to the /^-function coefficient in Feynman calculus in QCD. Thus, if one identifies 
the effective Hamiltonian form factor width parameter A with the running momentum scale in 
Feynman diagrams, the standard result from off-shell S-matrix calculus would be recovered in 
effective Hamiltonians in the case h{x) = 0. 

The function h{x) is determined by the initial Hamiltonian and small-x regularization factor 
rs in it through Eq. (4.13a). In the limit 6^0, for rg of Eq. (2.8), 



h{x) = 12 
and in the case of rs from Eq. (2.9) 



1 — X — x"^ , (1 — 3;)2 — a; , ,^ , 

3 + In X + — — In (1 - x) 

(l-a;)(l-2x) x{l-2x) ^ ' 



(5.6) 



h{x) = 12 In mm(a;, 1 — x). (5.7) 

These two cases are shown in Fig. 1 along with the case c) that corresponds to h{x) = 0. The 
vertex function V{x) in Fig. 1, is defined by 

goV{x)=go + Wx{x) , (5.8) 

and plotted using ao = go/{4:7r) = 0.1 for Aq = 100 GeV, and Nc = 3. In part d) of Fig. 1, the 
value of V{l/2) is plotted against A. The sharp cutoff case b) of 6{x — 6) is visibly different from 
c) that corresponds to h{xi) = 0, but the continuous case a) of x/{x + 6) differs from c) only 
by 8%. By extrapolation, the results a) and b) suggest that there may exist a regularization 
factor rs{x) that matches the case c). The mathematical structure of Eq. (4.13a) and Euclidean 
integrals in Feynman diagrams with dimensional regularization |^ both hint at the power- law 
functions rs{x) = x^. Inspection shows that in the limit 6^0, 
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a 



c 




b 




Fig. 1 Change of the effective gluon vertex function V{x) from Eq. (5.8) with the width 
A varying from Aq = 100 GeV down to 100 MeV, for three different small-x regularization 
functions: a) rs = x/{x + 6), b) rs = 6{x — 6) and c) rs = x^ 6{x — e). Part d) shows 
dependence of y(l/2) on A for the three cases, correspondingly. The case c, dashed hue, 
matches the QCD running coupling constant result obtained from Feynman diagrams. 
Note the dynamical suppression of the effective gluon coupling for extreme values of x in 
cases a) and b). See the text for details. 

rsix) = x^e{x-e) (5.9) 

leads to h{xi) = for all fixed values of Xi between and 1, as long as e/5 — > 0. This result sets 
a path for developing a connection between the Hamiltonian dynamics of effective gluons in the 
light-front Fock space, and Feynman diagrams for Green's functions in Lagrangian calculus [16]. 

The second scale e in Eq. (5.9) is required to regulate linear small-x divergences such as in 
the mass counterterms in Eq. (4.9). Details of rs{x) at x ~ e are not important for point-wise 
convergence of h{xi) to 0, and a whole class of regularizations with a second scale e <^ 6 give 
the same result. Ultraviolet regularizations using invariant mass, soften small-x divergences in 
mass counterterms so that the second scale e appears unnecessary, but such regularizations mix 
small-x effects with large k-^ divergences, see Appendix E. In evaluation of h{xi), factorization 
of the ultraviolet renormalization group flow in A from small-x effects requires S 0. The 
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small-x behavior of other terms in this hmit has to be compared with the size of g to discuss 
vahdity of the perturbative analysis. The warning is warranted by the fact that in Fig. 1, 
where the value of a ~ 0.1 for A on the order of 100 GeV is taken from phenomenology based 
on Feynman diagrams, ~ 1-1 > 1- 

Fig. 1 shows that small-x regularizations of Eq. (2.8) and (2.9) lead to suppression of 
effective gluon interactions when xi moves away from 1/2 for 0, where f\ 1. They 

also slow down the rate of growth of the coupling constant with lowering A. Similar results 
follow for fx depending on p~ instead of A4^. A^ is replaced by k^X and one looses boost 
invariance, but in the ratio (Xik^) / {X2k^) the dependence on ^3" cancels out. 

Thorn calculated the diverging part of a 4-point gluon Green's function using A'^ = 
gauge and integrating Feynman diagrams first over k~ . Using sharp cutoffs k'^ > e+ and 
\k^\ < A in the remaining integrals, Thorn found that to understand asymptotic freedom 
one has to include direct and crossed box diagrams that cancel skew p"*"- dependent terms 
from the three-point function whose own cutoff dependence had opposite sign to asymptotic 
freedom. Lepage and Brodsky |13[ developed a whole formalism for hard exclusive processes 
involving hadrons. Perry [^| used their rules in old fashioned perturbation theory to obtain 
the ultraviolet diverging part of o set of third order terms in the quark-gluon-quark off-energy- 
shell vertex function, with similar cutoffs to Thorn's. In the limit of ultraviolet cutoffs being 
sent to infinity. Perry found asymptotic freedom in the cutoff dependence of the quark-gluon 
vertex function. Irrelevant parts diverged as functions of the small-p"*" cutoff and depended 
on the off-shell energy parameter, but Perry reported that they could be made small by using 
an invariant mass cutoff. In calculations of the Green's functions, small energy denominators 
require an extra step of setting a lower bound on transverse momenta and factorization of the 
contributions from below. Such contributions complicate issues that arise in the introduction of 
the renormalization scale at small scales, and gluon coupling must be brought under control . 
The dependence on the energy-shell parameter involves Hamiltonian eigenvalues, which depend 
on solution of the hadronic binding problem. 

Effective Hamiltonians are calculated differently from Green's functions. Asymptotic free- 
dom appears in dependence of the Hamiltonian vertices on width A. The width may be finite 
and adjusted to the scale of processes one is interested in, independently of the diverging ul- 
traviolet cutoffs. Ux is unitary and the procedure involves no wave function renormalization, 
in distinction from the standard renormalization group concept, and no Ward identity is used. 
A effectively limits energy denominators from below, so that issues of binding are clearly sep- 
arated from logarithmic evolution of effective Hamiltonians. No off-shell energy parameter is 
introduced. In summary, no need arises to rely on the arbitrary ultraviolet regularization, no 
extra lower bound on momenta is needed, and no off-shell energy parameter is introduced. The 
notion of effective constituents is then a natural candidate for the phenomenology of hadronic 
wave functions |2i] to be put on the firm ground of Hamiltonian quantum mechanics, with an 
open path to make connection with diagrammatic techniques |^ for scattering amplitudes. 

Although the renormalized Hamiltonian vertices display a tendency to decouple dynamics 
of effective gluons from small-x region, many x dependent terms completely drop out from 
the third order running of the coupling constant. These terms require careful investigation 



since vacuum effects may enter through the small-x region |14]. Once these terms are cal- 
culated in detail using the perturbative procedure, they can be subsequently studied in the 
Schrodinger equations with effective Hamiltonians beyond perturbation theory. In particular, 
variational studies could aim at understanding gluon condensation and spontaneous chiral sym- 
metry breaking, including the cases where the number of flavors approaches the critical value 
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at which asymptotic freedom goes away |2B|. The small-x regularization clearly interferes in 
this transition, see Eq. (4.13). 



6. Conclusion 

Coefficients of products of creation and annihilation operators for effective gluons of width 
A in renormalized light-front QCD Hamiltonian, contain vertex form factors fx and vary with 
A in a perturbatively calculable way. Obtained in third order perturbation theory, the effective 
three-gluon vertex with vanishing transverse momentum, is a linear function of In A. The 
coefficient of In A matches the coefficient of In Q that appears in the running coupling constant 
dependence on the running scale Q in Lagrangian calculus for Green's functions. This happens 
for small-x regularizations of the type rs{x) = x^9{x — e) in the limit 5 — > for e/6 0. 
For other regularizations, such as rs{x) = x/{x + 6) or rs{x) = 9{x — 6), the coefficient of 
In A contains an additive term h{x), defined in Eq. (4.13a) and described in Section 5. h{x) 
suppresses interactions at small x in addition to the suppression implied by the vertex form 
factors fx. 

In the effective particle light-front Fock space basis in gauge A'^ = 0, asymptotic freedom is 
a consequence of that a single effective gluon contains a pair of bare gluons. This component 
amplifies the strength with which effective gluons can split into effective gluons when A gets 



smaller. The mechanism is the same as in the case of scalar particles in 6 dimensions [0, 
although the source of In A is different. Namely, in scalar theory it is the integration over 
additional transverse momentum dimensions, while in QCD it is the transverse momentum 
factors in gluon polarization vectors. In perturbative description of processes characterized by 
a physical momentum scale Q, using Hx{ax), there will appear powers of and InQ/A. For 
Q/\ = 1, \nQ/\ = 0, so that the theoretical predictions will have the form of a series in powers 
of the asymptotically free running coupling constant qq. 

The dynamics of effective particles includes binding, described through Schrodinger equa- 
tion with Hx{ax)- Initial studies of some simplified bound state eigenvalue problems for light- 



front QCD Hamiltonian matrices, have recently been carried out for heavy quarkonia |^ and 



gluonium states ||28[, reporting reasonable results. Although the effective particle approach 
described here is different, the bound state equations resulting from the second order perturba- 
tion theory may be similar. So, the initial studies suggest that the effective particle approach 
should be tested in application to bound states of constituent quarks and gluons. 

Besides the bound state dynamics of low energy constituents in the hadronic rest frame, 
the boost-invariant effective particle approach provides a theoretical tool for studies of layers 
of hadronic structure in other frames of reference, including the infinite momentum frame. The 
theory is simplest, and in lowest order of the same type as for quarks and gluons, in the case 
of electron-hadron interaction through one-photon exchange. In the case of massive vector 
bosons, there exist three instead of only two polarization states and the choice of gauge A'^ = 
is not directly available. A scheme to attempt in the more complicated theory was formulated 
by Soper in massive QED, but the required renormalization procedure for Hamiltonians 



with massive gauge bosons is completely undeveloped in comparison to the current status of 



renormalization in Lagrangian calculus |30|. 



Physical processes that involve one photon exchange have amplitudes proportional to e^. 
To obtain an amplitude for a simplest scattering event with large momentum transfer that 
involves strong interactions, it is sufficient to calculate the strong subprocess only to order g"^ 
and study cross-sections to order e'^g^. The complete initial Hamiltonian contains four terms. 



17 



H = Hqed + Heq^ + Hqcd + X , (6.1) 

where the second term includes couphngs of quarks to photons and instantaneous electromag- 
netic coupling of quarks and electrons. The transformation Ux of Eq. (3.1) can be calculated 
from the trajectory of Hxqcd alone, together with the required counterterms in X. Then, the 
Hamiltonians Hf,q^ can be re-written in terms of quarks and gluons of width A. This step pro- 
duces dependence on regularization, which should be removed by additional counterterms in X. 
QED degrees of freedom are not changed. The resulting effective Hamiltonian dynamics can 
be then tested for covariance using perturbative expansion for the amplitude e^e~ hadrons 
including terms order e^, e^g and e^g"^, that contribute to the cross section through orders 
and e^g"^. 

Such tests are of interest since the renormalization group for particles produces light-front 
Hamiltonians without reference to the vacuum structure and scattering amplitudes, and wave 
function renormalization and Ward identities do not enter the Hamiltonian. In a much simpli- 
fied model without gauge symmetry, adjustment of finite parts of counterterms led to covariant 
scattering with proper threshold behavior |^T|. It should be verified in the Hamiltonian ap- 



proach based on Eq. (6.1) if QCD binding effects known in order g"^ interfere with obtaining 
covariant results. Genuinely perturbative fourth-order test calculations in QED should verify 
if in the limit 5 — >^ the perturbative Hamiltonian approach is able to produce fully covariant 
results for observables through adjustment of finite parts of the ultraviolet counterterms that 
contain otherwise undetermined functions of x. 

Appendix A: Details of the initial Hamiltonian 

The initial QCD Hamiltonian for gluons in Eq. (2.11) contains the following terms. 

Ha^ =J2 J i^]-^^Lcakac ■ (Al) 
Ha^ = T. f [123] ^{p^ - p) rA5(3, 1) [g Yu3 aUla^ + g Y*^., ala2a,] , {A2) 

123 

where 



123 — 

with 6 = 6'^ and k = Hy^. 



1 . . 1 



^1^2 ■ ^3^^ ~ ^1^3 ■ ^2^3 ■ 

a;2/3 a;i/3 



{Ala) 



[1234] ^{p'^ ~ Y 1234444«4 + 1234aia2«3a4 + 12344^30201 



(A3) 



Sa4i234 = i[ ri+2,ir4,3 {elel ■ e\e^ - e\e^ ■ e*,el) f--^<^2 fac.c,^ i 

ri+3,1^4,2 {en ■ els, - els, ■ e*,el) /-i-3/ac,c4+ 3 / \4 ^ ^ 

^3+2,3^4,1 {en ■ ^2^4 - £3^4 ■ en) /-3C. jac,C4] , ^ 
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Xa^ 1234 = ri+2,ir3+4,3 (£^^3 " £3^4 - £^£4 " £^£3) ^^^^^ f<^-^-^ + 

[^3,1^2,4 + ri,3r4,2] (£^£2 • £3£4 - £t£4 • £2£3) /"'^i^3/ac2C4+ 1 — 

[r3,2ri,4 + r2,3r4,i] (£t£^ ■ £^£4 - £^£3 ■ £2£4) r^'^T^'^^ . 



3 
4 

1- 

2- 



-H"[aAA]2 = X] / [1234] (^(p^ -p) / (H[aAA]2i234ai44a4 + /i-c.) +X[aAA]2 1234440304 ■ (^4) 

1234 



^[BAAY 1234 = -\[ h+2,lh,3 SlS*, ■ SlS, ^"^^^f^^;"^^ fac-.c.fac,c,^ 
ri+3,1^4,2 £^£3 • £2£4 ^"iiZ^""'^ faccjacc,^ 
r3+2,3r4,l £^£^ • £^£4 ^"^^^jiy^^ /-ac^ jac,C4] . 

4[ n+2,ir3+4,3£l£2 ■ £3£4 (xi+xz)'^ ^ J ~ 

(xi+X;i)(X'2+Xi) 



[dAA]'^ 1234 



faC2C4 . 



[r3,if2,4 + h,3h,2] £^£3 • £^£4 ''''^::!}::p' r^^^r 

[h,2fl,, + f2,3f4,l] £^£4 • £^£4 ^"^fft;"^^ /aciC4yac2C3] . 



1 
2 
3 / 



-3 1- 

-4 2- 



and p are the total momenta of created and annihilated particles, respectively. For the vertex 
with parent momentum p and daughter momenta d and p — d, the regulating function in Eq. 
(A2) is fAs{p,d) = rAs{p,d)rAs{p,P - d), where r^sip^d) = rA{K,djl)r5{xd/p)0{xd/p) , with 
= d'^ /p'^ = Xd/xp and ac;^/^ = d"*" — XdjpP^ ■ Eqs. (2.6) to (2.9) complete the definitions of 
TAsip, d). Also, fp = rAsip, d). The momentum integration measures contain the same factors 
as in Eq. (2.4) for each indicated particle. 

Appendix B: Expression for 7;^2i 
Eq. (3.19) gives 

7A21 = ^3A 8 [Q;i2a2ia2l]21 + ^2\ 2 [/3oo22Cl;2l]21 + ^2\ 2 [Cl;i2/^oo3l]21 



+ Tz\ 4 [[q;i2Q;2i]iiQ;2i]21 + ^2\ [ai2/5oo3l]21 + ^2\ 2 [/3cx)lia2l] 
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+ ^3A 2 [a2l[Q;i2a2l] 11)21 + -^SA 2 [/^oo22a2l]2l + •^2A [a2l/3ooll]21 
+ 7oo21 ■ 

m 

The terms are grouped and ordered in one to one correspondence to Fig. 2, so that using labels 
from Fig. 2, Eq. (Bl) reads: 



721 = :S8(a)+JS2(6)+^22(c)+^34(d)+^2(e)+:S2(/)+^32(^)+:S2(^)+-^2(0 + (j) ■ (^2) 
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Fig. 2 Graphical illustration of Eq. (Bl). 

In all terms, the external gluon quantum numbers are labeled in the same way, 1 and 2 refer 
to the created gluons and 3 to the annihilated one. The intermediate momenta are assigned 
numbers 6, 7 and 8. For example, using notation from Appendix A, we have = Xipk^ + k^*^, 
where k^^ = K,j;^j, x^j-j = x^jx. Also, Kgg — — {\ — x^k^^I x-i^ k^q — —Xik-^/x + ^'^d /t^g = 
K"*". All coefficients of creation operators ought to be symmetrized with respect to quantum 
numbers of gluons 1 and 2, and only one way of contracting the intermediate gluon annihilation 
and creation operators is included, so that the displayed numerical weights represent the result 
of all possible contractions. Thus, in Fig. 2.b, only one ordering of vertices from /3oo22 is 
included and the factor 2 is put in front. The black dots in Fig. 2 indicate counterterms, /?ooii 
in cases f and i, and 7oo2i in case j. Thick lines with transversal bars denote the combined 
contributions of the terms H^a and H]^dAA]^ from Eqs. (A3) and (A4), /3oo22 = -^a* + ^[9aa]2 
and Poozi = + '^[dAA]^- Since these terms are independent of the transverse momenta and 
the three-gluon vertices ai2 and 0:21 are odd in the transverse momentum, both terms e and h 
in Eq. (B2), represented by diagrams e and h in Fig. 2, are equal zero. 

The non-zero terms are listed in the order of their appearance in Eq. (B2), using convention 
721 = En72i(n), with n ranging from a to j, and 

[123] ^{ki + k2- ks) - 7(„) alala^ . {B3) 
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7(a) = 

where 



dxrst{x) 



XI x{l — x)(x — Xi) 



±\ *^3A(a) i j k ijk , /-■ 



2) 



(S4) 



'^(5t(a^) = rs{x) rs{l - x) rs{xi/x) rs[{x - xi)/x] rs[{x - Xi)/x2] rs[{l - x)/x2] , (-B4a) 



TAti^^) = exp [-2(48^ + Kt,' + K^')/A'] , 



{BAh) 



3A(a) 



xMIq + M'^ 



fesfbdfief — 1 



{(Ml + x^Mls) 

f 68 fbd — 1 



Mis + ML + Mi, + M' Mi, + Mi 



'■bd 



'16 



68 



'■bd 



Mi, + M^ 



bd 



f 68 fief ~ 1 



f68fbdfl6f — 1 



Mi, + Mis . 

f fca 1 

' M"" + {M^ ~ Mj2)\ 
ffcafwfes 



Mis + Mi, + M^ Mis + Ml + Mi, + M^ 

ffcafiefes ~ 1 



(2M^ - M 



+ {M^ - Mi^f + Mi, + Mis 



Mi, 



Mis 



12 



//le/as — 1 



+ Mi, + Mi 



'68 



1 
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K 



±2 



M 



68 



x{l — x) ' 

2 I 2 
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{x — Xi){l — x) 



Ml, 



,.2, .±2 



(S4c) 
(BAd) 

(BAe) 
{BAf) 



X\{x — Xi) ' 

with = exp [-M^/A"^], hd = = Mls/x2 + M\2, ca = - MI2, f = fed, cd = M^, and 



Jjk _ *j *i k 
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El £2 ^3 



1 - 



X 



1 2x 
+ + 



X 



+ 



XX2 



+ 



X — Xi Xi Xi (1 — x)xi (1 — x)xi {x — Xi)Xi 



X — X\ \ — X 



I *k *j i 
-f- ti b2 £3 



-X2 
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X2 



XX2 



X — Xi {1 — x){x — Xi) 
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+£2^3 



iki ^2 



\l-xf \l-x)x 



k I XX2 X2_ 

^ ' — Xi)"^ 1 — X ^ 



X 



1 



XX2 



{1 — x){x — Xi) X ^ 



XX2 



(1 — x){x — XiY 



(1 — — Xi) 



*k 



{x — Xi)' 



X2 



x{x — Xi) 



XX2 



X2 



(1 — xYxi {x — Xi)x\ 



where 



2A(6) 



2M^ - MI2 



ifacf-l) , 



7(6) = ^ -y^^ y d K r^,(« ) ^ + (1 - 2) (S5) 



(B5a) 



and S{p) = (xi + x) (x2 + 1 - - xi)^. 



7(c) 
where 



^2A(c) _ X2Ml^+Mld 



•^2A(b) I n ^ 9^ 



^3^ 



A^^8 + Ml 



{fesfbd — 1) , 



(S6) 



(56a) 



^(c) = £(0)^^68 — ^1^2^3''^68 



-X2 X2S(c) 



X — X'l 1 — X 



and S(c) 



/ -, X2 X2 

+8^6382^68 -1 - ■S(c) + \ 

V 1 — X X — Xi 

{xi — X + xi){l — x + 



+ ^2^3 £1^68 



-X2 _^ X2S(^c) 



1 — X X — XiJ ^ 



{B6b) 



7(d) + 7(/) 



+41^123 



1 1 

1 + ^ + 



(1 — xy 



^3X{d) _^ ^2A(/) 



xlK X2k^M' 
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where 



{B7) 



3\{d) 
2.+ 2 



+ 



2A(/) 



1 



X2M' 



X2ktM^ xl Mi2 + M^ 

f fbdf 12 — 1 



M 



'i^-^^+X2j+M 

ffbd — 1 

2A^4 + (^Myx2 + Mh? + M\2 + {^1x2 + A^fs)'. 



12 



M'^/X2 



12 



/V: 



12 



1 



Pfbdf: 



12 



MI2 fl2p 



X2M^ M\2 + 2A^4 



(S7a) 

•^bd — M.^lx2 + AI12, since Algs = iii terms d and / (here / is a subscript, not a form 
factor), /12 = fad, ad = MI2, 



2A(/) 



/12-I 



^2^3" u.2^»'il2 

r^^Jyx) is given in Eq. (4.7), ri^^{K^) = rX{K^^), and 2g'^fij 



{B7b) 



167r3/i2. 



7(9) + 7(i) 



2iVcyi23 



Jo xil — x) J 



+2Y, 



123 



'0 x(l — x) 



1 + ^ + 



x^ 



{1-xy 



3A(s) , T2\(i) 



ktM^ 



{B8) 
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ktM 

f fca f \ 



12 



Mi2 + 2M^ M^ + Mi2 

- 1 ffca - 1 



{(2M'-A< 



2 ^ 
12 J 



2M^ + (M^ 
+ M\2 



^\2, 



12 
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M2 



ff 



12 



1 



f fcaf\ 



A^2_^2 



42J 



12 



2M^ + Mi2 2M^ + {M^-Ml2Y + M\2_ 



with ca = M'^ — MI2, and 



2X{i) 



f. 



12 



kt 



12 



7(j) = 7oo + (1 <-> 2), where 700 denotes the counterterm coefficient in 



{B8a) 



(586) 
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[123] ^{ki + k2- ka) 700 a\ala3 , {B9) 

whose dependence on the quantum numbers 1, 2 and 3 needs to be found. 

Appendix C: Ultraviolet divergent terms in 7^21 

The three-gluon vertex is given by Eqs. (B4-9). The uhraviolet diverging parts of terms 
b and c vanish. This can be seen by changing variables so that the sum of squares of relative 
transverse momenta in three subsequent vertices in Eq. (B4b), is written as 

V6s48+Vi6i^w+Vi^^^^CP^^ + Xi^i2^ , (CI) 

where C = Ves+Vw^l/^'^+V: = «"^-C«i^2> C = [V68{'i--3c)/x2+rji6Xi/x]/C: and x = [^68^i6(a^- 
/ {xx2y + r)Qsv{^ ~ 1^2 + ^i6^]/C- The coefficients rj arc introduced for identification of 
finite parts of counterterms that may contain functions of Xi. All 77s are equal 2 in Eq. (B4b). 
Using variable p^, the potentially diverging relative transverse momentum integrals in terms h 
and c, can be written as 

Hp)<uv = -2 ^^^3 1^ ^^^^^(^) f ^2^± g-Cp-/A^ Jl_ + (1^2) , (C2a) 

li,c)aiv = 2 ^i/--c3 I' ^^^^^(^) I ^2p± g-Cp-/A^ + (1^2) . {C2h) 

Z Jx\ J Kgg X2 

Both terms contain an integral of the same structure, 

= j dt'K^ ^ exp [-C{k^ - v^f/^^] , (C2c) 

and limA^oo I'^ — t^v^, which is a finite ultraviolet regularization effect with no divergence. 

Thanks to the regularization factors rst, these terms can be integrated over x and the resulting 
function of Xi depends on the choice of the coefficients 1]. If they were chosen to depend on x, an 
arbitrarily complex, ultraviolet regularization dependent finite function of Xi can be obtained, 

but ^(b)div = l(c)div = 0. 



^ {g+i)div 

-NcYi2z [ dxrs^i{x) 
Jo 



2 K,^ 



g-4«VA^ X(l - X) 



1 + ^ + 



x^ (1 — x) 



+ (1^2) . 



(C3) 

yU,^ is an arbitrary parameter inserted here to simplify notation for the well defined integral over 
at small k when the terms with form factors are written explicitly and the small denominator 
effect is absent. One also obtains 



l{d+f)div — '^l[g+i)dv 
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The logarithmically diverging term (C3) appears then with factor —3 in the countcrtcrm co- 
efficient 7oo in Eq. (B9). Its finite contribution to the counterterm is a finite number times 

^123- 

The remaining term, 7(0) of Eq. (B4), contains the tensor k*-'*^ = kI^KiqK^ contracted with 
£^^^ and the latter does not depend on transverse momenta. The ultraviolet regularization 
dependence comes only from 



•^3A(a) 



X2 



and using Eq. (CI) one obtains 

1{a)div ~ 



-^ir^'^--^ / dxutix) +(1^2) , 



where 



pjk ^ J ^2 pi. g- 



AC ACgS 



(C5) 



(C6) 



(C7) 



The ultraviolet-regularization-dependent part of P^'^ can be evaluated using Feynman's trick 
to combine denominators and making similar steps as in Eq. (C2), 



pjk 



71 Xil — X 

^ 2 X X2 



i $:jk k ~^ XX2 j r-n^ 

t\j'\f)0 t\j'\t)0 — 1~ ^ ^12 



— x) 



In 



l'^12l 



(C8) 



A number of finite terms are obtained that depend on the choice of coefficients r] in Eq. (CI). 



The terms include also a new tensor structure Ki2'^i2'^i2/'^i2- Inserting 
in Eq. (C6), one obtains 



pjk 



in place of 



ijk 



where 



and 



with 



-"C.- fClC2C3 



l{a)div^^'-^ir"'''' dxrstix) 

Z Jxi X2 



1 — X 



pjk 



+ (1 - 2) 



a"(«) 



if 



C1C2C3 



1 — X 

X2 



C12 



jijk 



6?. =71 In — [C12Y12 + VisYis + C23Y23] 



2 1 1 (l-x) 
+ + -+ ^ ^ 



1 — X X — Xi X 



X2 



!_ A 

X2 



Cl3 



2 1 1 (1-xf ^ 
+ + - + ^ '--2 



1 — X X — Xi X 



X2 



2 1 

C23 = \- 



1 — X X — Xi X 



1 (1-xf l+x^ 
+ --^ ^ 2 , 



x^ 



X2 



(C9) 

(C9a) 

(C96) 
(C9c) 
{CM) 



Y,2 = ip'^'^elel ■ £3^12 , 



(ClOa) 
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Y23 = -ir'^'^'ele^ ■ £t«i2 



1 

2^2/3 
1 



The counterterm should cancel the divergence, so that 

= l{a)div + l{d+f)div + l{g+i)div + loodiv + (1^2) 



A r fl 

2Nc'K In - <^ 2 / rfx r^t (x) [012^12 + CisYis + 0231^23] 



-3 [yi2 + ^13 + 1^23] C dxrs^ix) 
Jo 



x(l — x) 



1 



1 + .2 + (1 _ ^)2 



and 



where 



'y{j)div = -^123^^^ In — [11 + h{xi)] , 

O jJi 



{cm) 

(ClOc) 



+ loodiv + {I ^'2) , (Cll) 



(C12) 



/i(a;i) = 6 / dxrst{x) 



2 11 

H h 



-9 / dxrsJx) 
Jo 



1 — X X — Xi X 
The ultraviolet counterterm ^0021 has the form (B9), where 



1 1 

+ 



X 1 — X. 



+ (1^2) 



Too = ^123—^ In — [11 + h{xi)] + Ifinite ■ 
o 1^ 



(C13) 



(CU) 



Different choices of the coefficients rj in Eq. (CI), lead to different finite terms 1 finite: which 
demonstrates finite dependence on ultraviolet regularization. To restore covariance of observ- 

ables, the finite terms must be then allowed to contain unknown functions of Xi that multiply 
three structures Y12, F13, and Y23 from Eq. (CIO), and the fourth structure. 



(C15) 

Inclusion of n/ flavors of quarks produces in 'j(j)div additional diverging terms of the form 



-27m f In — [cfuYu + 0/13^13 + C/23I23 + CfmYua] , 
A* 



where 



c/12 



C/13 



1 1 - 4a; + 2a;2 + xi 

7^ + (1^2) , 

XI ^X2 



— dx 

Jx\ 



1 , 1 - 4x + 2a:^ + xi 



2X2 



+ (1^2) , 



/•I -2 + 4x-2x2(l + a;2) +a;ia;2 , ,^ ^, 
c/23 = - j^^ dx ^ + (1^2) 



(C16) 

{Clla) 

{cm) 

{Cllc) 
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rl _3 

Cfm = / dx—2[l - 2x{l - x)] + (1^2) . (C17d) 
^0 8 

After simplifications, integration and symmterization, fermions' contribution turns out to be 
not sensitive to small-x regularization and changes Eq. (C14) to 

lij)div = -^23^ In - {iVe [11 + h{x,)] - 2nf} . (C18) 



Appendix D: W\{x) in Eq. (5.2) 

W^Ai2(a^i, '^12^), Wxi3{xi, Kji2^) and Wx23{xi, in Eq. (4.14) become all equal to Wx{xi) 
when Kj;2 = 0. Calculation of W\{xi) is based on extraction of coefficients of terms linear in 
, which form Y^j for ij = 12, 13, 23, and employs the following facts. The self-interaction 



K, 



12) 



terms d and g, and mass counterterms / and i from Fig. 2, contribute through 



167r3 



Nr / dx 







d^ 



K 



-_L2 



■x{l — x) 



1 + ^ + 

x^ 



1 



[i-xy 



X 



m 

II is introduced only for simplification and is canceled in the full formula for dependence of 
Wx{xi) on A. In the term 7(a), the renormalization group factor of Eq. (B4c) contributes to 
the dependence of Wx{xi) on A only through small A expansion of 



3A(a) 



= -[Mh,-2x,A]M^^^''^'' ' ^^^^ 

where A — K^K^2/ix — Xi) and the added subscript indicates that is set equal 0. Other 
parts do not contribute because when the form factors / are expanded, no dependence on A is 
generated due to dimensional reasons and the remaining factors cannot contribute since they 
are multiplied by differences of form factors and there is an identity 



i 



°° dz 



2 a 



which shows that A drops out. The expansion in A in Eq. (D2), produces the same tensor 
structure as in Eq. (C8), which leads then to Eq. (5.2). 



Appendix E: Regularization mixing for x and k 

Regularization with invariant masses implies that the coefficients ri in Eq. (CI) become 
equal r\ = 1/x + 1/(1 — x), r^es = X2/{x — xi) + 0:2/(1 — x) and rji^ = x/xi + x/{x — xi). 
These coefficients diverge when x ^ xi or x ^ 1. In mass counterterms, the coefficient 77 is 
the same. When integrating over transverse momenta in mass counterterms with rAi2'r)K-^^), 
one obtains in Eq. (4.9) only logarithmically divergent integral over x and a regularization 
factor of the form rs{x) = x^ is sufficient to regulate it. In other words, in place of the second 
scale, e, in rs from Eq. (5.9), one can consider regularizations of transverse divergences that 
provide additional damping of the small-x region. The question is what are the consequences 
of the mixing of large k, and small x in ^X2i- Most representative is Eq. (B4). Although 
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all three rjs can grow to oo, x is limited and does not exceed l/{xiX2) reaching this value at 
the ends of the integration range, while ^ = 1 when x = Xi and drops down to at a; = 1. 
Therefore, x^\2 I ^ i^i the exponent is always a small number and vanishes when A — cx) 
without contributing to regularization dependence. On the other hand, the coefficient C grows 
to infinity at the ends of the integration region in x. This way the ultraviolet regularization 
factor depending on invariant masses changes small-x singularities. The same phenomenon 
occurs in a simpler form in the instantaneous terms involving 1/9^^. The instantaneous terms 
themselves do not contain integration over k^, but the integration occurs when these terms are 
included in the dynamics. The logarithmically divergent integrals over produce logarithms of 
C as A-independent remnants of the ultraviolet regularization and one cannot exclude arbitrary 
functions of x in finite parts of the ultraviolet counterterms, including integrals that strengthen 
small-x logarithmic divergences. Invariant mass regularizations including a small gluon mass 
11 P], provide additional damping for x ~ yU^/A^. 

References 

[1] Particle Data Group, Eur. Phys. J. C3, 1 (1998). 

[2] R. P. Feynman, Phys. Rev. Lett. 23, 1415 (1969); Photon- Hadron Interactions (Benjamin, 
New York, 1972). 

[3] See e.g. M. A. Shifman, A. I. Veinstein, V. I. Zakharov, Nucl. Phys. B147, 385, 448 (1979). 

[4] K. G. Wilson, Phys. Rev. DIO, 2445 (1974). 

[5] Nucl. Phys. (Proc. Suppl.) B83-84, Lattice '99 (2000). 

[6] St. D. Glazek, Acta Phys. Pol. B29, 1979 (1998); Phys. Rev. D60, 105030 (1999); Acta 
Phys. Pol. B31, 909 (2000). 

[7] St. D. Glazek and K. G. Wilson, Phys. Rev. D48, 5863 (1993); zbid. D49, 4214 (1994). 

[8] See e.g. K. Erkelenz, Phys. Rep. 13C, 5 (1974), Sec. 2.3; St. Glazek, Acta Phys. Pol. B15, 
889 (1984). 

[9] K. G. Wilson et al., Phys. Rev. D49, 6720 (1994). 

[10] J. B. Kogut and D. E. Soper, Phys. Rev. Dl, 2901 (1970); J. B. Bjorken, J. B. Kogut and 
D. E. Soper, ibid.3, 1382 (1971); and references therein. 

[11] S.-J. Chang, R. G. Root and T.-M. Yan, Phys. Rev. D7, 1133 (1973); S.-J. Chang and 
T.-M. Yan, ibid, bf 7, 1147 (1973); T.-M. Yan, ibtd.7, 1760 (1973); ibid.7, 1780 (1973); 
and references therein. 

[12] S. J. Brodsky, R. Roskies and R. Suaya, Phys. Rev. D8, 4574 (1973). 

[13] G. P. Lepage, S. J. Brodsky, Phys. Rev. D22, 2157 (1980); and references in there. 



28 



[14] J. Kogut and L. Susskind, Phys. Rep. C8, 75 (1973) and references therein. For a more 
recent presentation, see L. Susskind and M. Burkardt, in Theory of Hadrons and Light- 
Front QCD, edited by St. D. Glazek (World Scientific, Singapore, 1995), |hep-ph/94103l3|; 
K. G. Wilson and M. Brisudova, zbid., |hep-th/9411008| ; St. D. Glazek, Phys. Rev. D38, 
3277 (1988). 

[15] St. D. Glazek and K. G. Wilson, Phys. Rev. D57, 3558 (1998); cf. F. Wegner, Ann. Physik 
3, 77 (1994). 

[16] D. J. Gross and F. Wilczek, Phys. Rev. Lett. 30, 1343 (1973); H. D. Politzer, ibtd., 1346. 

[17] H. J. Melosh, Phys. Rev. D9, 1095 (1974); and references therein to a series of earlier 
papers that contributed to the formulation of QCD. 

[18] Second paper in 0. 

[19] G. Altarelh, G. Parisi, Nucl. Phys. B126, 298 (1977). 

[20] See e.g. M. E. Peskin, D. V. Schroeder, An Introduction to Quantum Field Theory 
(Addison- Wesley, 1995), Eq. (7.80). 

[21] C. B. Thorn, Phys. Rev. D20, 1934 (1979). 

[22] R. J. Perry, Phys. Lett. B300, 8 (1993). 

[23] S. J. Brodsky, G. P. Lepage, P. M. Mackenzie, Phys. Rev. D28, 228 (1983). 

[24] S. J. Brodsky, H. G. Pauh, S. S. Pinsky, Phys. Rep. 301, 299 (1998); S. J. Brodsky, 
|hep-ph/ 00092281 and 229. 

[25] A. Bassetto, G. Nardelli, R. Soldati, Yang-Mills theories in algebraic non-covariant gauges 
(World Scientific, Singapore, 1991). 

[26] T. Appelquist, Z. Duan, F. Sannino, Phys. Rev. D61, 125009 (2000); and references 
therein. 

[27] M. Brisudova and R.J. Perry, Phys. Rev. D54, 1831 (1996); M. Brisudova, R.J. Perry and 
K. G. Wilson, Phys. Rev. Lett. 78, 1227 (1997). See also R. J. Perry and K. G. Wilson, 
Nucl. Phys. B403, 587 (1993) and R. J. Perry, Ann. Phys. 232, 116 (1994). 

[28] B. H. Allen, R. J. Perry, Phys. Rev. D62, 025005 (2000). 

[29] D. E. Soper, Phys. Rev. D4, 1620 (1971). 

[30] G. 't Hooft, M. Veltman, Nucl. Phys. B44, 189 (1972); ibid. B50, 318 (1972). For reviews 
of further development, see G. 't Hooft, Nucl. Phys. B (Proc. Suppl.) 74, 413 (1999); D. 
J. Gross, ibid., 426. 

[31] T. Maslowski, M. Wi§ckowski, Phys. Rev. D57, 4976 (1998). 



29 



